High frequency vortex dynamics and magnetoconductivity of high temperature 

superconductors 



Jan Kolacek 

Institute of Physics, ASCR, Cukrovarnickd 10,16200 Prague 6, Czech Republic 

Etsuo Kawate 

National Research Laboratory of Metrology, 1-1-4, Umezono, Tsukuba, Ibaraki 305, Japan 

(February 1, 2008) 

The vortex lattice with the superconducting and normal state charge carriers fractions may be 
regarded as three independent subsystems mutually connected by interactions. The equations of 
motion for these three subsystems must be solved simultaneously. In this way a new consistent 
theory of the vortex dynamics is obtained and the magnetoconductivity calculated. 
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I. INTRODUCTION 



The vortex dynamics attracts great attention from both the theoretical and experimental point of view. Much 
controversy is concentrated on the question, whether Magnus force is the only transverse force on the vortex as is 
claimed by Ao and Thouless or if other ones, like Iordanskii force from the normal fluid or Kopnin-Kraftsov force 
from impurities may also contribute to the total transverse force. This topic is discussed in detail by Sonin ||, where 
many relevant references can be found. 

As the characteristic frequencies of the vortex system are in the far infrared (FIR) region || the magnetooptical 
spectroscopy is a suitable tool to study these problems. Recently results of the FIR magnetotransmission of YBaCuO 
thin films using elliptically polarized light became available pH. To interpret the measurements, the high frequency 
conductivity must be calculated. The first such calculation was made by Gittleman-Rosenblum |(|. Later, the 
dissipation near the flux line depinning threshold using the generalized temperature assisted flux flow (TAFF) model 
was studied by Yeh . In these papers the reaction of the superconducting (and normal) charge carriers on the vortex 
movement was neglected. The redistribution of the induced ac current density due to the vortex response is taken 
into account in the Coffey-Clem |E] calculation of the surface impedance. Other references concerning this problem 
may be found in a recent review [ 3 1 . 

In principle, our approach is similar to the one of Coffey-Clem model, but the formalism is different. We treat the 
superconducting fluid, normal fluid and the vortex system as three subsystems mutually connected by interaction. 
Their equations of motion are solved simultaneously. In this way a new, internally consistent theory of vortex motion 
is obtained. We use this approach for calculating magnetoconductivity in high temperature superconductor (HTSC) 
thin films. In this paper we consider the Magnus force as the only interaction between the vortex system and 
superconducting fluid and Lorentz force to be the main interaction force between the normal fluid and vortex system 
in HTSC materials. It conforms to the necessary symmetry requirements summarized by Sonin [Q] and has the same 
sign and form as the D,D' terms derived by Stone ]io| for normal particles confined in the vortex core, if the limit of 
nonlocalized particles is taken. We show, that including this transversal force from normal state charge carriers gives 
reasonable magnetoconductivity tensor which is compared with recent experimental data [Q . 

II. THE MODEL 



Let us consider the vortex lattice, in which the mean distance between vortices is small in comparison with the 
penetration depth. In this case the magnetic field in the superconductor is almost homogeneous. This model may 
be appropriate for HTSC. In the field 5 T the mean distance between vortices is about 20 nm, while in YBCO the 
London penetration depth Al is of the order 150 nm. The diameter of the vortex core is approximately equal to 
the coherence length £, which in YBCO is about 2 nm. Therefore, it is possible to neglect possible redistribution of 
the charge and current density in the vortex core and to take it into account only as the source of vortex damping. 
The damping force is supposed to be proportional to the vortex velocity with frequency independent viscosity 
coefficient r\, so that 
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F D = -VVl = v L . (1) 

T„ 

Introducing m v as the vortex mass per unit length, we will use the vortex relaxation rate 1/t v which is given in 
practical frequency units, rather than the viscosity. 

For pinning we use the simplest model of a parabolic well, so that the pinning force Fp is proportional to vortex 
displacement r^: 

Fp = -nr L = -m v a 2 r L (2) 

with k and a being the pinning constant and pinning frequency, respectively. 

The interaction between superconducting fluid moving with velocity v s and the vortex system is mediated by the 
Magnus force jjj given by 

ft h 

F M (v) = -y-(v a - v L ) x z = m v f s n(v s - v L ) x z, (3) 

where n s = f s n is the density of superconducting fluid and Fjf(n) means, that this is the force felt by the vortex. 
The reaction force Fjvf(s) acting on a superconducting particle is 

F M (s) = -F M {v) = -mu c (-v s - Vi) x z, (4) 

n s 

where n v is the vortex density (number of vortices per unit area), u c = eB/m = n v h/2m is the cyclotron frequency 
in the field B = n v &oz caused by the vortex system and $o is the magnetic flux quantum. In (|^) we introduced the 
frequency of the cyclotron vortex motion Q = nh/2m v . It is interesting to note that, using the Hsu's expression for 
the vortex mass m v = (7r 2 /4)mfc F £ 2 H], the 2D expression for the Fermi wave vector k 2 F — 2itn and £ = Tivf/kIS. 
for the coherence length, it is possible to show, that fl = A 2 /Ep (A is the gap and Ep is the Fermi energy), which 
is the level separation in the vortex core M. 

The interaction between the vortex system and the normal state fluid may be obtained in the following way. From 
the Aharonov-Casher Lagrangian (l2) it can be shown that if the vortex lattice moves with velocity vl, the force 
imposed by the vortex system on one normal state particle moving with velocity v„ is 

7i h 

F L(n) = -^-( v » ~ v l) x z = mw c (v B - v L ) x z. (5) 

According to the action-reaction law the vortex lattice must feel the same force with opposite direction. If there are 
n n = f n n normal state particles, the total force per unit length of one vortex is 

F L (v) = -— Fi(n) = -/ri^r(v„ - v L ) x z = -m„/ n fi(v„ - y L ) x z. (6) 

This expression is analogous to the Magnus force formula, but has opposite sign. It satisfies the invariance require- 
ments, according to which only the relative velocity of the particle with respect to the vortex system is decisive. Factor 
/„ is justified by the fact that the total force is proportional to the number of particles involved in the interaction 

The questions concerning forces acting on the vortex lattice are still seriously controversial. The Lorentz force 
(^) following from the Aharonov-Casher Lagrangian is of electrodynamic origin, but similar formula is also used to 
describe interaction of normal state fluid with vortices in neutral systems (see e.g. [§10 i1|) . Useful comments and 



replies regarding the spectral flow force and the Iordanskii force can be also found in |14-|16|. Interaction of electric 
charge with moving vortex and the Aharonov-Casher effect in two-dimensional superconductors was discussed e.g. by 
Simanek fPfl ] . Let us note that it would not be correct to consider the Lorentz force also for the superconducting fluid, 
which would exactly cancel the Magnus force. The vortex lattice and the accompanying magnetic field are created by 
superconducting current, so in this case Lorentz force would mean "action on itself. 

Having draw up the interaction forces we will now write the equations of motion for the three subsystems. As in 
the London model, the superconducting fluid is supposed to move without damping, 

mv, = eE + F M (s), (7) 
while the normal state fluid motion is damped as in the conventional Drudc model 

mv„ = eE + F L (n) v„. (8) 

Tn 
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Finally, for the vortex system we shall use the Newton type equation of motion (of course the vortex mass and also 
all the forces are considered per unit length) 



m v v L = F P + F D + F M (v)+F L (v). (9) 



The system of three equations of motion (|7||9j) together with expressions for the interaction forces (|3|-|6|), damping 
and pinning force (QJ|) form a closed set of equations for the unknown v^, v s and v n . Assuming a periodic time 
dependence e luJt , the three differential equations reduce to the set of three linear equations : 

A ss v s + A sv vl = — E 
to 

A nn v n + A nv v L = —E (10) 
to 

v s + A v n + A vl = 

with the coefficients 

A ss = i(u - lo c ) ; A sv = iuj c 

A nn = i(w + uj c - i/r n ) ; A nv = ~iuj c 

A vs — if s Q J A vn = —ifrSl (11) 

A vv = i(u> + (f n - f a )Cl - a 2 /u - i/r v ). 

If the determinant D = A ss A nn A vv — A nn A sv A vs — A 8S A nv A vn is nonzero, the set can be readily solved to get 

gE A n nA vv -\- A sv A V n A n vA vn cE 



v s = g s 



m Dm 



cE A SS A VV -\- A n vA vs A SV A VS cE . 

V n =9n = Fi (12) 

TO U TO 

cE A ss A vn A n nA vs cE 

v L = g L = . 

to Dm 

Now it is straightforward to express the conductivity as 

cr = = ^(n s v s + n n v n ) = e ^l(f s g s + f n 9n), (13) 

where iv p — \J ne 2 /e^m is the plasma frequency and the factors g s , g n are defined by eq. ([l2|). 

As expected, for physically meaningful parameters (r v > 0, r„ > 0, Lu c tt > 0) the real part of conductivity 
is positive and the Kramers-Kronig relation ct(uj) = ct + (w/i7r) j_ 00 <j{x)/(x 2 — xu>)dx as well as the f-sum rule 

(I/71-) J_ oo Re(a(cj))cLo — €quj 2 are satisfied. The zero frequency limit of the conductivity Co = eoUpT n (f n +i[u) c T n (f s — 
f n ) + fs/uj c T n ])/(l + t 2 u> 2 ) does not have the delta function component, as the pinning constant is supposed to be 
finite, while pinning range is infinite. Expressing the conductivity tensor components as a xx (u>) — (c(u)) +<j(— u>))/2 , 
o xy (u}) = (cr(ui) — a(—u>))/2, it is possible to show, that also the Hall sum rule pl| (1/tt) f^^Re^H )duj — u>h, where 
tn = & X y/&xx and loh = lim w ^ 00 [— iLutn (to)] = uj c (f s — /„) is satisfied. It is necessary to note, that without normal 
state fraction ( /„ = 0) the tn function has a pole at zero frequency, so that in this case the Hall sum rule must be 
modified to lor = (a 2 LU c /(QuJc + a 2 )) + (l/V) J^ c oo Re(tH)du;. 



III. ABSENCE OF NORMAL STATE FLUID 



Usually it is considered that at zero temperature all the charge carriers condense, so that normal state fluid is 
absent. It is not necessary true for all materials, but it is useful to discuss this limit first. 

For free vortices (vortices without pinning and damping) the two equations of motion tov s = — tow c (v s — vt) x z 
for the superconducting fluid and tovl = m v f s Q(v s — vj,) x z for the vortex system are readily simplified to vl/v s — 
Q/(p, — bj) and vl/v s — (lu c — lo)/oj c , respectively. Consequently, two nontrivial solutions exist: for zero frequency 
Vl = v s , while for u = Q + oj c the velocity ratio is vl/v s = —Q/lo c . This means that the superconducting liquid and 
vortices may move either in parallel with constant velocity (this solution is required by Galilean invariance), or may 
oscillate with opposite phase, with the inertial center remaining at rest. 

In general, with /„ = the coefficient A vn equals zero and the conductivity formula ([l3]) reduces to 



3 



a(f n = 0) = e ^ A ™ . (14) 

Let us note that in the limit of zero vortex density (lo c — > 0) , this formula reduces to the London expression for con- 
ductivity a = eoWp/iuj, as expected. For zero pinning (but nonzero damping) the explicit expression for conductivity 
may be written as: 

*(/» = 0, a = 0) = VT ^? V ( 15 ) 

r„w(w c + si — cj) + i(w — lo c ) 

It is clear that in this case the real part of conductivity is nonzero even at zero frequency <Ji{f n — 0, a = 0, lu = 0) = 
eou>pT v n/uj c . Contrary to it, for nonzero pinning we get a(f n — 0, a ^ 0, lu = 0) = eoLu 2 i/Lu c with zero real part of 
conductivity. This result is understandable, if we recall that in our simple model the pinning barrier is infinite, so 
that the d.c. transport must be nondissipative. 

In reality the pinning barrier is not infinite. Depending on frequency, temperature, magnetic field, as well as density 
and strength of pinning sites, various regimes as flux creep, flux flow, temperature assisted flux flow etc. fl9[| can be 
recognized. To keep the discussion simple, we will analyze just two simple limits. In the "full pinning" (FP) limit 
the driving field is low, so that each vortex is bound to the individual pinning valley, making only small oscillations. 
In this case the pinning force plays an important role. Contrary to it in the limit of high driving field the amplitude 
of the vortex oscillation is larger than the distance between the pinning centers, and the averaged pinning force is 
effectively zero (ZP). In the intermediate state the pining force is nonzero, but not proportional to the distance from 
the pinning center which leads to nonlinear effects. We will show that, in some frequency range, nonlinear effects can 
be expected even at relatively low fields which are commonly used in laboratory experiments. 

Let us estimate the realistic values for the parameters of the theory. For coherence length £ = 2 nm, and effective 
mass m = 4m e using the Hsu's expression for the vortex mass jy], we can estimate il = 2h/(n 2 mt; 2 ) — 49 cm -1 . 
The cyclotron frequency in the field 4T is 5.9 cm -1 . Using the expression || k = (0.01 -r- 0.05)iiqH 2 for the pinning 
coefficient, and the vortex mass estimation |t| m v = 1.6*10 10 m e /m, the range for pinning frequency a = 19-^95 cm -1 
may be obtained. As we did not select any model for the vortex damping, we leave this parameter as free. To make 
a model calculation we used the following set of parameters: lu c — 5, £1 = 50, a — 30, lu p = 6000, l/r v = 10 (all 
values are in cm -1 ). The conductivity calculated in FP and ZP limits are displayed in fig.l. The conductivity peaks 
are expected near the frequencies, where the real or imaginary part of the determinant D — A SS A VV — A SV A VS which 
appears in the denominator of ( |l4| ) is zero. In ZP limit the expected peak values of conductivity are 

<ti(u = 0) = e Q ujpT v fl/uj c 
a\ {to = fl + w c ) = eoUjpT v u; c /n (16) 
<ti(w = uj c ) = e uip/T v Vtuj c . 

In fig.l. only two sharp peaks are present for the ZP limit (dashed line). It is obvious from ( |l6|) that, while for 
low vortex damping (large t v ) the peaks at eigenfrequencies of the system (0 and Q + uu c ) are important, for high 
vortex damping the peak at lu c will dominate. This is illustrated in fig. 2, where the conductivity for vortex damping 
1/t„ from 10 to 200 are displayed. It is possible to see, how with increasing vortex damping the peak shifts from zero 
frequency to the cyclotron frequency uj c . For FP limit, due to the pinning term a 2 /lo the order of the determinant D 
is higher in u>, so one more peak is expected in accord with the model calculation results displayed in fig.l (solid line). 

In fig. 3 the relative value of the vortex oscillation amplitude a v = \vi\mjeE is shown as a function of lu. It is clear 
that, while at high frequency the oscillation amplitude is low so that FP limit is appropriate, at lower frequencies the 
amplitude is high, so that ZP limit must be adopted. In principle, beside the pinning frequency a determining the 
pinning force at low oscillation amplitude, two characteristic lengths r\ and V2, the amplitudes of vortex oscillation, 
at which the pinning force declines from the linear law and at which the pinning force is effectively zero, should be 
introduced. In this way, for a given driving field E, two crossover frequencies Wdi,d 2 = \JeE /n^m are defined. For 
a v < the FP limit is valid, while for a v > l/^ 2 the ZP must be used. If neither condition is fulfilled, the 

system is in a nonlinear region, where the conductivity depends on the driving field strength. It is interesting to 
note, that for some frequencies both conditions a v (FP) < l/w^i an d av(ZP) > l/io^ 2 may be fulfilled at one time. 
This means, that in this frequency region bistability may occur. Depending on the history, at the same experimental 
conditions two regimes - the low and high vortex oscillation amplitude, corresponding to the low and high resistivity 
state - may be achieved! All these possibilities are illustrated in fig. 3. If we estimate the range of pinning force (r^) 
to be about 10 nm and if the intensity of radiation used for measurement is 1 mW/mm 2 so that the driving field E 
is of order 1.7 * 10 4 V/m, we get LOd = 9 cm -1 . For illustration purposes we have chosen = LOd 2 = 10 cm -1 . It is 
obvious that, depending on the frequency and intensity of the radiation used for the measurements, many interesting 
nonlinear effects may be expected. 
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IV. INFLUENCE OF NORMAL STATE FLUID 



For nonzero temperatures there are two contributions to the real conductivity. One is connected with the normal 
state charge carriers, the other with vortices. As expected, without vortices we get <t(uj c = 0) = eaUp[f s /iw+f n T n / (1 + 
iujT n )\ , which is the sum of the London and Drude model contributions. The normal state limit (f s — > 0) does not 
have much sense, as without superconducting fraction we can not have any vortices. However, if we simulate external 
magnetic field by making vortices unable to move, we should get the formula for a normal conductor in magnetic field. 
Indeed, in the limit of vortices fixed to the lattice (a — > oc) or of infinite vortex mass (f2 = 0), we get the expected 
result a(f s = 0, O = 0) = a(f s = 0, a — > oc) = e uj^T n / (1 + i(u + w c )r„). 

The results of model calculations for /„ = 0.5 in FP and ZP limits are displayed in fig. 4. To visualize the contribution 
of vortices, the zero magnetic field conductivity (u> c = 0) is also displayed in these graphs. In FP limit, when the 
amplitude of vortex motion is small, almost all real part of conductivity originates from the normal state charge 
carriers - except of the very sharp feature near the zero frequency, which is caused by the vortex resonance. On the 
other hand, in the ZP limit the conductivity is much larger and it is almost completely due to the vortex motion, 
with the normal state fluid playing only a minor role. However, the sharp vortex resonance peak is absent. It might 
be somewhat surprising, that the presence of vortices may slightly decrease the real part of conductivity for some 
frequencies. 

It is instructive to see, how increasing the normal state fraction influences the conductivity. For the FP limit it 
is shown in fig. 5. We can see that with increasing /„, the central peak (connected with the normal state carriers 
conductivity) gradually develops, while the side peaks diminish. Recently, Lihn at al. {| measured far infrared 
magnetoconductivity tensor in YBaCuO thin film. Their data are also displayed in fig. 5 (dashed line). The intensity of 
FIR radiation is usually rather small so the FP limit could be appropriate. It is remarkable, that all the experimentally 
observed features are quite well simulated by the curve with /„ = 0.3. This seems to indicate the presence of some 
normal state fraction (probably located on CuO chains ) even at the lowest temperature. Alternatively it may be due 
to the enhanced density of quasiparticles in an applied magnetic field, as predicted for the d-wave superconductor 
|^,|lj. It should be noted, that the sharp vortex resonance peak on the theoretical curve is at lower frequency than 
the range accessible by FIR spectroscopy, so it could not be observed in the experiment. 



V. CONCLUSIONS 



Vortex lattice together with the superconducting and normal state fluid form three subsystems mutually connected 
by interaction. Taking into account reaction forces by which vortices influence superconducting and normal state fluid 
and solving simultaneously the three equations of motion a new, internally consistent theory of vortex dynamics was 
developed. It was shown that due to the finite range of the pinning force, at some frequencies nonlinear phenomena 
may be expected even for relatively low driving fields which are commonly used in laboratory experiments. For 
comparison with experiment, the knowledge of the power of radiation used for the measurements might be crucial. 
The presented theory can qualitatively explain recent measurements of far infrared magnetoconductivity tensor made 
by Lihn at.al The d.c. conductivity calculated in the framework of this model enables to explain theoretically 
controversial, but experimentally firmly established Hall voltage sign reversal 
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FIG. 1. Real part of conductivity for FP (solid line) and ZP (dashed) limits. Parameters of the model are lo c — 5, Q = 50, 
a = 30, lo p = 6000, \/t v = 10 (all in cm" 1 ). 



FIG. 2. Real part of conductivity for ZP limit, for uj c = 5, fi = 50, a = 0, uj p = 6000 and 1/t„ = 10, 20, 50, 100, 200 (all in 
cm -1 ). For increasing vortex damping the conductivity peak shifts from zero frequency to cyclotron frequency. 

FIG. 3. The relative vortex oscillation amplitude for FP (solid line) and ZP (dashed) limits calculated using the same 
parameters as in Fig.l. The pinning range corresponding to crossover frequency uj d = 10 cm -1 is marked by a horizontal line. 
The frequency regions marked by \\\ (///) are regions where FP (ZP) limit are appropriate. 

FIG. 4. The comparison of real part of conductivity in the FP (solid line) and ZP (dashed) limits with conductivity in zero 
magnetic field (dotted) in presence of normal state fraction /„ = 0.5; l/r„ = 15. All the other parameters are same as in the 
Fig.l. 

FIG. 5. Real part of conductivity in FP limit for normal state fraction /„ = 0.1, 0.15, 0.2, 0.25, 0.3 (the tendencies for 
increasing /„ are marked by arrows). Other parameters are same as in Fig.l. The experimental data obtained by Lihn at al. 
[4] are also plotted (dashed line) for comparison. 
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